extern doublereal dlapy3_(doublereal *x, doublereal *y, doublereal *z__);
